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Abstract. 

We compare the two-year COBE DMR sky maps with the predictions of cosmological- 
constant cold dark matter (ACDM) models. Using a Bayesian analysis, we find that the 
most likely value of the cosmological constant in a ACDM model is A = 0. The data set 
an upper limit on A of 0.78 (0.85) at 90% confidence, and 0.86 (0.92) at 95% confidence 
with (without) the quadrupole anisotropy. The angular power spectrum C; depends only 
quadratically on A when A is small; the likelihood function L(A) is therefore quite fiat 
near A = 0. 

Subject headings: cosmic microwave background — large-scale structure of universe 
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Introduction. 

Over the last few years, both the quantity and the quahty of cosniological data have 
improved rapidly. Several new observations and data analyses of the large-scale structure 
of the Universe have been performed (see e.g., Peacock & Dodds 1994). Moreover, the 
recent detection of anisotropy in the cosmic microwave background (CMB) by the COBE 
satellite (Smoot et al. 1992) provides important information about the initial conditions 
of the fluctuations. It is now possible to test models for the growth of large-scale structure 
in the Universe more precisely than ever before. The general picture that appears most 
consistent with these data is that structure formed by gravitational instability from small 
perturbations in the early Universe. Within this paradigm, however, there are many 
unanswered questions. The manner in which structure evolves depends upon the amount 
and kind of matter in the Universe, as well as the statistical properties of the initial 
fluctuations. The most popular model so far is the cold dark matter (CDM) model, in 
which the Universe is assumed to be cosmologically flat, with density parameter il = 1. A 
small fraction (around a few percent) of the matter in this model is baryonic, while the rest 
is nonrelativistic coUisionless matter of some kind. The initial perturbations are assumed 
to be Gaussian and adiabatic, with a power-law power spectrum P{k) = k"^, with n sa 1. 

The CDM model's popularity arises from the fact that it is approximately consis- 
tent with observations of galaxy clustering, and with observations of the CMB anisotropy. 
However, as the quality of these data has improved, including the detection of the CMB 
anisotropy by the COBE DMR (Smoot et al. 1992), signiflcant discrepancies have arisen 
between the observed power spectrum of galaxy clustering and the CDM prediction (Efs- 
tathiou. Bond, & White 1992, Peacock & Dodds 1994). In order to resolve this problem, 
there have been a variety of proposed modiflcations to the CDM model: among the pro- 
posed flxes are tilting the primordial power spectrum (see e.g., Davis et al. 1992), mixing 
some amount of hot dark matter with the cold (see e.g., Davis, Summers & Schlegel 1992) 
and assuming a small value for the Hubble constant (Bartlett et al. 1994). One very 
promising modiflcation involves replacing some of the cold dark matter with a cosniologi- 
cal constant A while keeping a flat geometry (see e.g., Efstathiou, Sutherland & Maddox 
1990; Efstathiou, Bond & White 1992; Kofman, Gnedin & BahcaU 1993). This has the 
added advantage of solving the cosmic age problem. 

It is known, however, that A models produce quite different CMB anisotropies on 
large scales (Kofman & Starobinski 1985), because of the time evolution of the gravita- 
tional potential which produces the Sachs- Wolfe integral (Sachs & Wolfe 1967). Sugiyama 
& Silk (1994) numerically solved the perturbation equations and obtained CMB power 
spectra for different A's. The power spectra for A-dominated models have shapes that are 
quite different from the simple power-law models that are usually used in analyzing the 
CMB anisotropy. Therefore, these power-law analyses (e.g., Smoot et al. 1992, Seljak 
& Bertschinger 1993, Gorski et al. 1994) cannot be used to set limits on A. A direct 



comparison between the DMR data and the predicted power spectra from the A models is 
required. 

In this paper we compare predicted CMB power spectrum of ACDM models with the 
COBE DMR data. In these models, the Universe is supposed to be cosmologically flat, 
so that the density parameter il and the cosmological constant A add up to unity. The 
matter in the Universe is assumed to consist primarily of CDM, with a few percent of 
the matter being baryonic. It should be noticed that the CMB anisotropy on the angular 
scales probed by COBE is largely independent of the matter content of the Universe. 
Therefore our results obtained here are general for ACDM models. The initial fluctuations 
are supposed to be Gaussian and adiabatic, with a Harrison- Zel'dovich (n = 1) power 
spectrum, as in standard CDM. These models flt the galaxy clustering data quite well 
with A ~ 0.8 (Efstathiou, Sutherland, & Maddox 1990; Efstathiou, Bond, & White 1992; 
Kofman, Gnedin & Bahcall 1993). 

Method. 

We compare the DMR data with the ACDM predictions using a likelihood technique 
similar to that described by Bond (1994a, 1994b). We will assume that the CMB anisotropy 
seen by the DMR has the properties of a Gaussian random fleld. It is useful to write the 
anisotropy as an expansion in spherical harmonics: 

AT v-^ 

-jri^) = Z^(iimYiUr). (1) 

l,m 

(Here and hereafter, all sums over / run from 2 to oo, and all sums over m run from — / to /. 
Modes with / = 0, 1 are marginalized over during the likelihood analysis.) The coefficients 
aim are independent Gaussian random variables with zero mean and variances 

C, = (|a,^|2), (2) 

where the angle brackets denote an ensemble average. The power spectrum C; is related 
to the rms temperature fluctuation in the following way: 

(|AT/Tp> = 5](2/ + l)C,/47r. (3) 

The angular power spectrum C; is determined by the cosmological model under consid- 
eration. Our goal is to determine which models (i.e., which angular power spectra) are 
consistent with the DMR data. 

The DMR data consist of Np = 6144 pixels. After excising all pixels within 20° of the 
Galactic plane, the number of pixels remaining is N = 4038 (for maps that are pixelized in 
ecliptic coordinates). Each pixel contains a measurement of AT/T at a point r^, smoothed 



with the DMR beam pattern, and contaminated with noise: the zth data point can be 

represented as 

/AT \ 
di = ( — ;- * B ] (ri) + ni = y^ BiaimYimiri) + rii, (4) 



l,m 

where B represents the beam pattern and the star denotes a convolution. Bi is the Zth 
coefficient of the expansion of the beam pattern in Legendre polynomials. It can be 
approximately modeled as a Gaussian, Bi sa exp(— 0.5/(/ + 1)/19.3^), but we use the 
more accurate values from Wright eA al. (1994). Ui is the noise in the zth pixel. To a good 
approximation, the noise is uncorrelated from pixel to pixel (Lineweaver eA al. 1994), and 
the quantities Ui can be accurately modeled by Gaussian random variables with known 
variances cr?. 

Our goal is to assess the consistency of various cosmological models with the data di. 
Ideally, we would like to compute the likelihood of each model, 

-2lni: = (FM-^(f+lndetM, (5) 

where the covariance matrix M is determined by the angular power spectrum C; predicted 
by the model: 

Mij = {didj) = J2 BfCiYUri)Yrjrj) + a]Sij. (6) 

l,m 

Unfortunately, it is impractical to compute L for each model we wish to test: 4038 X 
4038 matrices are not easy to invert. We are therefore led to consider ways of approximating 
the likelihood function. 

Since the DMR beam pattern drops off rapidly with /, we expect that a relatively 
small number of modes, say those with / < Imax, contribute significantly to the spherical 
harmonic expansion of AT/T. It therefore seems plausible that we could "compress" the 
data vector in some judicious way to a dimension D ~ /^^j^, that would be significantly 
smaller than N . Then the likelihoods would be much easier to compute. If we performed 
this compression in the right way, we could hope that the compressed data would contain 
most of the useful information in the raw data. If the data covered the sky completely, 
then the obvious way to effect this compression would be to use the data to estimate the 
individual spherical harmonic coefficients: bim = {4:TT /N)Y^-_-^ diYim{i^i)- In the case of 
incomplete sky coverage, the orthogonality of the spherical harmonics is destroyed, and 
this procedure no longer provides reliable estimates of the true coefficients (Peebles 1980; 
Bunn, Hoffman & Silk 1994). However it is still perfectly acceptable to use these quantities 
for computing likelihoods, as long as one takes proper account of their covariance matrix 
(Gorski 1994). 

We will follow a slightly different approach. In order to avoid repeated inversion of 
large matrices, we will compress the data in the following way before computing likelihoods. 

— * 

We will choose some D X N matrix A with D < N , and compute u = Ad. Then for any 
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model we wish to test, we can compute the D X D covariance matrix of u, M^ = AMA, 
and the hkehhood, 

-2lni:„ = M^M-^M + lndetM„. (7) 

In order for these hkehhoods to be useful in constraining models, we need to make 
a good choice for the matrix A. In particular, we would like to choose A in such a way 
that the likelihood function L„ has maximum rejection power. Suppose the correct power 
spectrum is C;. Then if we make a small change to the power spectrum, say by changing 
the normalization by a factor (1 + e), we would like the likelihood to decrease as much as 
possible on average. That is, we would like to minimize (SLu) = {Lu{Ci) — Lu{{l + e)C;)). 
Since (Lu) has a maximum at the true model e = 0, (SLu) oc (L'l), the second derivative 
of (Lu) with respect to e evaluated at the true model e = 0. 

It is straightforward to calculate (L'l), and to minimize this quantity with respect to 
the matrix A. The optimal matrix A turns out to be a matrix whose rows are solutions a 
of the generalized eigenvalue problem 

M,iga = \Ma, (8) 

where the covariance matrix M consists of a signal part and a noise part: Mgig^ij = 
Y^i,m^l^l^lrn{ii)Yi*^{ij) and Mnoise,ij = (^f^ij. Furthermore, the eigenvectors we want 
are those with the largest eigenvalues A. 

This result is not difficult to interpret. The rows of A are, to use Bond's terminology 
(Bond 1994a, 1994b), "eigenmodes of the signal-to-noise ratio." The components of the 
reduced data vector u are simply the projections of the raw data vector d in the direc- 
tions of these eigenmodes. The modes with large eigenvalues are the directions in the 
iV- dimensional data space which have large sensitivity to the signal, as compared with the 
noise. In this sense, the signal-to-noise eigenmodes are the optimal directions to look at 
in data space. 

Of course, in order to determine the eigenmodes, we need to make a guess about 
the actual power spectrum of the data. However, it turns out that there is little danger 
in guessing wrong: the likelihood function L„ does not change significantly as the input 
power spectrum is varied. For definiteness, from now on we will take the rows of A to be 
the solutions of (8) with a Harrison- Zel'dovich power spectrum. 

We also need to decide on a value for D, the dimension of the reduced data vector 
u. We want D to be large enough to have good rejection power, but small enough to be 
computationally convenient. We have found by trial and error that the likelihood functions 
for D = 400 and D = 500 are almost identical, whereas the D = 300 likelihood function 
is noticeably broader, signifying less rejection power. We have therefore chosen to use 
D = 400 in testing our models. With D = 400, it is possible to compute all rows of A 
by truncating all sums over / at Imax = 30: increasing Imax to 40 leaves these eigenmodes 
unchanged. 



Although we remove a best-fit nionopole and dipole from the data, in the absence of 
complete sky coverage we cannot be sure we have removed precisely the correct monopole 
and dipole. We therefore marginalize over these quantities by integrating over all possible 
values of the monopole and dipole. This is the correct procedure when performing a 
Bayesian analysis; an alternative procedure would be to take the peak likelihood as the 
monopole and dipole are allowed to range over all possible values; this procedure gives 
results which are negligibly different from marginalizing. 



Results. 

We have applied the technique presented in the previous section to the COBE DMR 
two-year sky maps. The DMR data consist of six sky maps, two independent channels at 
each of three frequencies. We formed one sky map by taking a linear combination of all six 
maps, with weights proportional to the inverse square of the noise levels in the maps. We 
excised all pixels within 20° of the Galactic plane, and removed a best-fit monopole and 
dipole from the data. We then computed the signal-to-noise eigenmode matrix A and the 
reduced data vector u for this data set. It is now straightforward to compute the likelihood 
Lu for any model we wish to test. 

Before examining ACDM models, we tested our method on standard power-law power 
spectra. In an il = 1 cosmology, on angular scales large enough that the only significant 
source of anisotropy is the Sachs- Wolfe effect (Sachs & Wolfe 1967), the angular power 
spectrum is 

„ f4.Q^\ r(ii±pi)r(V) 



5 J r (^^±1^) r (^) ' 

where the power spectrum of density perturbations is P{k) oc k^ and the quadrupole 
Q = Qrms-PS = Y 5C2/47r is used to normalize the power spectrum (Abbott & Schaefer 
1986; Bond & Efstathiou 1987). We computed the likelihood function Lu{n,Q) and used 
it to constrain the slope n and normalization Q of the power spectrum. We found that 
the maximum-likelihood solution was (n, Q) = (1.3, 17.1 //K). Furthermore, using standard 
Bayesian techniques with a uniform prior, we found that 0.64 < n < 1.78 at 95% confidence, 
and that 17.9 < Q < 24.3 for n = 1 at 95% confidence. Both results are consistent with 
previous analyses (Gorski eA al. 1994). 

Figure 1 shows the angular power spectra C; for seven ACDM models. These power 
spectra are obtained by direct numerical calculations up to the present epoch (see Sugiyama 
& Gouda 1992). All have Hubble constant iJo = 50kms~-^ Mpc"-^ and baryon density 
Q,B = 0.03. The cosmological constant A takes the values 0,0.2,0.4,0.6,0.7,0.8,0.9. (The 
model with A = is, of course, standard CDM.) The dependence on Q,b of C; is quite 
weak, i.e., less than 10% for / ^ 30. The dependence on the Hubble constant is also weak: 
changing Hq from 50 to 80 leaves the C;'s almost identical for / ^ 10, and causes less than 



a 20% difference for / ^ 30. It is clear from Figure 1 that the C;'s for these models are 
quite different from flat or simple power-law models, particularly when A is large. For each 
of these models, we computed L„ as a function of the power spectrum normalization Q. 
These likelihoods are shown in Figure 2. It is clear that the data prefer low values of A. 

We can place Bayesian confldence limits on A in the following way. For each of the 
seven models, we compute the marginal likelihood by integrating over the normalization 
Q: 

L^arg{A)= I L4A,Q)dQ. (10) 

These marginal likelihoods are plotted in Figure 3 as a function of A. We then smoothly 
interpolate between the data points, using a natural cubic spline. We can then say that A 
has an upper bound Amax at some confldence level c if 

/ L^arg{A-)dA = C L^arg{^)dA. (11) 

Jo Jo 

The upper limit on A derived in this way is, to use the standard statistical nomenclature, 
the boundary of a Bayesian credible region in the Q-A plane, adopting a prior distribution 
that is uniform in both Q and A (Berger 1985). A uniform distribution is the proper "non- 
informative prior" for A. Several different choices of prior distribution could be justifled 
for Q- however, since the data contain a very strong detection, the results are insensitive 
to the prior in Q (Bunn et al. 1994). Table 1 shows the upper limits for A for several 
different confldence levels. 

An alternative procedure would be to use the maximum value of the likelihood, 
Lmax{^) = niaxQL„(A, (5), in place of the marginal likelihood in eq. (11)- In prac- 
tice, it makes very little difference which procedure is followed: replacing L^arg by L^ax 
changes the 95% confldence level upper limit on A from 0.86 to 0.87. 

The quadrupole moment of the DMR sky maps is anomalously low, and it is possible 
that it may be contaminated in some way (see e.g., Gorski eA al. 1994). We explored the 
sensitivity of our results to this possibility by removing the quadrupole from the data and 
repeating the analysis including only modes with / > 3, marginalizing over the quadrupole 
as well as the monopole and dipole. Table 1 contains the upper limits on A found in 
this way. Since the large quadrupole predicted by ACDM models is at odds with the low 
quadrupole in the data, it is not surprising that removing the quadrupole weakens the 
constraints. 

The likelihood curve in Figure 3 is flat near A = 0. The reason for this is easily seen 
by looking at Figure 1: the angular power spectrum for A = 0.2 is almost identical to the 
A = model, indicating that C; is a very weak function of A for small A. Speciflcally, if 
we perform a Taylor expansion of C; about A = 0, the linear term is negligibly small. 

Conclusions. 
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We find that ACDM models are strongly disfavored by the DMR data. In particular, 
models with A > 0.86 are inconsistent with the data at a confidence level of 95%, while 
models with A > 0.78 are inconsistent at the 90% level. The large-scale structure data 
seem to prefer ACDM models with A ^ 0.7. It is clear from these results that the DMR 
data are at best marginally consistent with such models. 

Removing the quadrupole from the data weakens these constraints somewhat. How- 
ever, it is important to note that the primary reason for removing the quadrupole is simply 
that it is anomalously low. Removing data points simply because they do not fit our theo- 
retical expectations is a dangerous statistical practice. In the absence of strong independent 
evidence that the quadrupole is contaminated, one should be wary of throwing it away. 

These results we obtained here are about the same as the best current limits from 
gravitational lensing (Kochanek 1993, 1994). It should be noticed that our method is 
totally independent and free from any ambiguities, such as the evolution and mass distri- 
bution of galaxies, involved in the gravitational lensing method. There is every reason to 
expect that, once the full four-year DMR data have been analyzed, it will be possible to 
discriminate between models with A sa 0.6 and those with A = with very high confidence, 
since the power spectra in these models have very different shapes (see Figure 1). 
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Figure Captions. 

1. The angular power spectrum /(/ + 1)C; is shown for each of the seven ACDM models 
described in the text. The power spectra have been normalized so that 6C2 = 1. 

2. The likelihood L„ is plotted as a function of Q for each of the seven ACDM models 
described in the text. The overall normalization is arbitrary. 

3. The marginal likelihood L^arg is plotted as a function of A. The solid curve is the 
likelihood including the quadrupole, and the dashed curve is the result of excluding the 
quadrupole. The crosses show the marginal likelihood for each of the seven models de- 
scribed in the text. The smooth curves are cubic spline interpolations between them. The 
overall normalization is arbitrary. 
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